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ABSTRACT 


An important goal of a tidal theory is the improvement of nutational amplitude and of the 
parameters of the Earth’s elastic response. With this goal in mind, a theory of tidal oscilla- 
tions inside a rotating elliptical Earth is developed, with special emphasis on tides in the 
liquid core. The Molodensky and Kramer theory of the resonance effect, as caused by the 
proximity of the frequency of the free diurnal wobble of the liquid core to the frequency 
of Kj astronomical tide, was amended to include the effect of the possible deviation of the 
liquid core from the state of neutral stability. Coupling effects between the toroidal and 
spheroidal oscillations, as caused by the Coriolis force, are taken into consideration. 
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BASIC NOTATIONS 
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Stress in the Earth as induced by the lunisolar tidal forces 
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Total stress 
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Constant of gravitation 

Initial force function of self-gravitation (per unit of mass) at 
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ON THE TIDAL OSCILLATIONS OF THE 
LIQUID CORE OF THE EARTH 


Peter Musen 

Goddard Space Flight Center 
Greenbelt, Maryland 


INTRODUCTION 

From a parametrized model of the Earth, as obtained from seismic information, a model of 
the tides in the interior of the Earth can be determined. The understanding of the internal 
tidal processes is important from the astronomical and geophysical standpoints because they 
influence the amplitudes of nutation and the parameters of the Earth’s elastic response (Love 
numbers), as observed on the Earth’s surface. Thus, an important goal of a tidal theory is the 
improvement of the nutational amplitudes and of the Earth’s elastic response— pure, uncon- 
taminated by the effects of the free diurnal wobble of the liquid core. With this goal in mind, 
a theory of tidal oscillations inside a rotating elliptical Earth was developed, with emphasis 
on tides in the liquid core. It is planned to apply this theory to the actual computation of 
the tides, particularly in the outer core. 

The present work can be considered as an extension (or a modification) of tidal theories by 
Molodensky (References 1 and 2), Molodensky and Kramer (Reference 3), and Shen and 
Mansinha (Reference 4). 

Some of the final equations are identical to the final equations of Shen and Mansinha, and 
some are not. The differential equations for tidal amplitudes were deduced. Some of the 
equations by Shen and Mansinha clearly represent linear combinations of the equations for 
amplitudes. 

In the exposition the vectorial and dyadic symbolism is being used, partly to obtain the 
differential equations in a compact form and partly to clarify the kinematics of numerous 
couplings between the tides of different degrees, as well as between the spheroidal and 
toroidal tides. The basic equations include the effect of the Earth’s deviation from the 
hydrostatic equilibrium in the initial state and the influence of the vorticity on the stress 
tensor, making use of Biot’s (Reference 5) and Dahlen’s (Reference 6) formulations of the 
elasticity theory. 

Satellite observations indicate the existence of a small deviatory prestress. Biot’s and 
Dahlen’s works pave the way for its inclusion when it is considered opportune. 



The tides and tidal forces were expanded into a series of complex vectorial harmonics, both 
spheroidal and toroidal (References 7 and 8). This expansion presents no numerical compli- 
cations because modern electronic machines are now able to handle functions of complex 
arguments without difficulty. Some relationships were also established between vectorial 
harmonics of different degrees. These relationships facilitate the expansion of the Coriolis 
force into a series of vectorial harmonics and, consequently, the formulations of the differ- 
ential equations for tides. Further, they eliminate the customary use of 3j-Wigner symbols, 
which are not convenient for this problem. As in the work of Shen and Mansinha (Reference 
4), the original Molodensky and Kramer equations (Reference 3) were amended by including 
the perturbative effects of the ellipticity of equipotential surfaces and the deviation of the 
liquid core from the neutral stability. 

The established equations indicate precisely when and why the terms that contain the 
ellipticity as a factor can be neglected. Although the influence of the ellipticity on the 
astronomical tides is relatively small, it affects the free diurnal wobble of the liquid core. 

It must also be taken into account in the transfer of the normal component from the liquid 
core to the mantle through the ellispoidal core-mantle interface. The period of free diurnal 
oscillations of the Uquid core is only 3 minutes short of the period of the Earth’s rotation 
(Reference 4). This proximity of periods causes a strong resonance of the free diurnal oscilla- 
tions of the liquid core with the diurnal astronomical tide Kj and, to a lesser degree, with 
the diurnal astronomical tide, Pj (References 9 and 10). The influence of this strong reso- 
nance contaminates the observed diurnal Love numbers. The removal of this contamination 
is an important geophysical problem because it helps to obtain the Earth’s elastic response 
pure to the astronomical tidal forces only. 

The resonance effect is one possible cause of the dependence of diurnal Love numbers on 
frequency. Observations with horizontal pendulums and gravimeters lead to different 
values of Love numbers for different diurnal tidal constituents (Reference 10). Recent 
values of Love numbers for the “whole Earth,” as deduced from laser observations of 
GEOS-3, also display a marked dependence on frequency (References 1 1 and 12). Haardeng- 
Pedersen’s computations (Reference 13), based on the model of a rotating Earth with the 
effects of nutation and diurnal wobble of the liquid core excluded, also show the dependence 
of Love and Shida numbers on frequency. 

Poincare found the existence of the free diurnal wobble of the elliptical liquid core (Reference 
14), and Jeffreys discussed the geophysical implications of it (Reference 15). In their theory 
of the wobble, Molodensky and Kramer (Reference 3) assumed the neutral stability of the 
outer core and the validity of the Adams-WiUiamson condition (Reference 16). Shen and 
Mansinha, Haardeng-Pedersen, Pekeris and Accad, and Crossley (References 4, 13, 17, and 
1 8) recently studied the dynamical consequences of the deviation of the liquid core from 
the neutral stability. 

In the present work, Molodensky and Kramer’s theory of resonance (Reference 3) was 
amended by including the perturbative effects proportional to the index of stability of the 
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liquid core. Perturbations in the geopotential and gravity are controlled by the Poisson 
equation. If the Earth is spherically symmetric, then only spheroidal tidal components 
enter into the Poisson equation and produce a change in gravity and, because of the absence 
of lateral inhomogeneities and dilatational changes, toroidal oscillations do not perturb the 
gravity. In the elliptic liquid core, however, there are lateral variations in density, and 
toroidal amplitudes enter into the Poisson equation with the ellipticity as a factor. Thus, 
the toroidal oscillations can also produce the perturbations in geopotential and in gravity 
if they are sufficiently large and if there are marked lateral inhomogeneities in the structure 
of the Earth. 

Computation of tides in the mantle and in the liquid core are essentially two different prob- 
lems. In the mantle, tides are basically static, and only spheroidal tidal oscillations are easily 
observable. By comparison, the toroidal tidal oscillations in the mantle are small and are 
usually neglected. In the liquid core, the rigidity vanishes or is very small. As a result, the 
rotation of the Earth (Reference 19) and the ellipticity of equipotential surfaces induce 
toroidal tidal oscillations, as well as numerous couplings. 

The differential equations given in this work contain the effects of the Coriolis force, of 
the ellipticity of equipotential surfaces, and of the couplings between the toroidal and 
.spheroidal tidal oscillations, as well as the influence of the possible departure of the outer 
core from the neutral stability. 

It is of interest to determine the perturbative influence of viscosity of the liquid core on 
tides and on amplitudes of nutation. However, this influence cannot be easily estimated. 
Information on the viscosity in the outer core is incomplete, and the estimates of the kine- 
matic viscosity coefficient, v, vary in a large interval. At the present time, the building of 
tidal models for different values of v can only be attempted. 

RECURSIVE FORMULAS FOR ASSOCIATED LEGENDRE FUNCTIONS AND 
THEIR DERIVATIVES 

This section contains some recursive relations between the associated Legendre functions 
and their derivatives. These relations are used in performing the expansion of tidal oscilla- 
tions and tidal forces into a series in vectorial harmonics. 

Beginning with 


xP^(x) = 


n - m + 1 
2n + 1 




n + m 
2n + 1 


pm 

n-1 


(X) 


( 1 ) 
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the following is derived: 


and 


From 


. (n - m + l)(n - m + 2) . . 

x^P”(x)= .. . . 


(2n + 1) (2n + 3) 

2n^ - 2m^ + - 2n - 1 
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and, taking equation 1 into account. 


(1 - 


dx 


n(n - m + 1) 
2n 


c, « 


(» t l)(n ■> .n) . 
2n + 1 


Equations 1, 2, and 5 yield: 
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x(l - x^) " - 

dx 
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From the normalization conditions, 

+ 1 

pm(x)pm(^)clx = 

where 



2 (n + m)! 

” 2n + 1 (n - m)! 

and 6^^ are Kronecker’s deltas, and, taking equation 1 into consideration, 



X P|^ (x) (x) dx = + 
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2n + 1 
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and taking into account 
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Multiplying equation 12 by x and integrating by parts yields: 


li 


(1 . x^) “CW . <1CW ^ 
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— [n (n + 1) + p(v + 1) - 2 ] ^ X P™ (x) PJJ* (x) dx 


and, taking equation 10 into account, 
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(1 - x") 


2. dP^(x) 
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From equation 12, taking equation 8 into account, 
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2 ^ dP™(x) 
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From equations 6 and 8, it follows that 
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(1 - x^) 


dP“(x) 


P“ (X) = + 
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(n + 2) (n + m + 1) 


2n + 3 


“n ^.n+1 


(n - l)(n - m) 
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oc 8 - 
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SPHERICAL VECTORIAL HARMONICS 

In this work, the tidal oscillations and the tidal forces are expanded into a series in vectorial 
harmonics (References 3, 7, 8, 20, and 21): 


A«(0,¥>) = = e^e^i^-^P^Cx) 


(15) 


where 


B|;’(5,v’) = rVS^ 


[ dP*" txi ~\ 
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(17) 
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X = COS^ 


(18) 


= i sin 0 cos (/? + j sin 0 sin + k cos 9 


e = + i cos 0 cos (/? + j cos 0 sin (^ - k sin 6 


= - i sin V? + j cos ip 
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(19) 


and 


*r ^ ^ " V V ^ *r " ®e 


d Id 1 d 

^ 3r r 90 r sin 0 difi 


( 20 ) 


In the frame of the present (linearized) theory, there is no mutual influence between vectorial 
harmonics of different orders. For this reason, the index m in the exposition will be omitted 
in all cases in Mdiich ambiguity does not arise. 


From equations 15 through 17, taking into account 

= - n (n + 1) r'^ 


( 21 ) 


r . VV S = - V S„ 

n n 
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yields 
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(33) 


Making use of 
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the following is derived from equations 1 5 through 1 7 : 
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The asterisks designate complex conjugates. The use of equations 38 through 46 facilitates 
the expansion of the Coriolis force into a series of vectorial harmonics. In performing the 
expansions of tides and tidal forces into a series of vectorial harmonics, the inner product 
is used: 


(F, G) 



(47) 


where F and G must be of the same order. The dot designates the standard scalar product. 
From equations 13b and 15 through 17, the following normalization and orthogonality 
conditions are derived: 


(A„> A*) = % 5„„. 

^ (B„,C„.) = (C„,A„.) ^ 0 


(48) 

(49) 

(50) 


and from equations 38 through 46, 8, 10, and 12 through 14, the following relations are 
derived: 


(k X A„, A„,) = 0 

(51) 

(k X A„, B„,) = - ima 
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(55) 
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(kX B„,C„.) = - 
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(59) 


which are useful in the process of expanding the geostrophic force. 

The expansion of a vector, F, in terms of spherical vectorial harmonics has the form: 


(60) 


where all vectors are of the same order. On several occasions, the series on the right-hand 
side is finite, and it will be simply a finite linear combination of vectorial harmonics. In 
particular, taking equations 5 1 through 59 into account. 


im 
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DIFFERENTIAL EQUATIONS OF TIDAL OSCILLATIONS 

Let V be the volume of a portion of the Earth in the “initial” state. Assume that in v and on 
its boundary surface, a, the density and elastic parameters are continuous functions of the 
position vector, r. 

Let dv be the element of v, and let da be the oriented surface element of a. At the moment, 
t, under the influence of tidal forces, these quantities become V, A, P, R, dV, and dA, respec- 
tively. Assume that the tidal displacement 

u = R - r 

is small and that all quantities of the second order in u can be neglected. From 


1 

dA = — dR A dR 
2 


9R 

dR 

dR 

dR 

dR 

dR 

X 

by 

dz 

dy dz + 
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X dz dx + 

dx 

dx 

X — 
dy 


( 64 ) 


and 


da 


£ 

2 


dr A dr = i dy dz + j dz dx + k dx dy 


( 65 ) 
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the following is derived: 


dA = da • A 


where the dyadic A has the form: 


3R 3R 9R 3R 3R 9R 

A = i X + i X + k X 

9y 9z 9z 9x 9x 9y 


Substituting 


R = r + u 


into equation 67 and taking the identities, 

IX i = kj-jk, IXj=ik-ki, IX k = ji-ij 

into account, then neglecting the terms of the second order in u. 


Taking 


A = I - (I X V) X u 


(I X V) X u = uV - IV • u 


and 


uV=e + IX fa> 


( 66 ) 


(67) 


( 68 ) 


into consideration, 


A = (l+d)I-e-IX<o (69) 

Under the influence of vorticity, the prestress, S, will receive at the moment, t, the increment 
(Reference 5) 

CO X S - S X CO 


and the total-stress tensor becomes: 

T = t + S + co.X S-SX« 


(70) 
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where T is the stress induced by lunisolar tidal forces. Let U (R) be the force function of 
self-gravitation (with centrifugal effects excluded), and let W be the tidal-force function. Let 
(j> be the change in the gravitational-force function as caused by a small redistribution of 
matter in the Earth’s interior under the influence of tidal forces. 

In this case, the principle of D’Alembert takes the form: 




V [U(R) + i//] - w|PdV + / dA*T = 0 


/ 


(71) 


where 


= W + 0 


a^R 3R dJ2 

w = + 2i2 X + X R X (n X R) (72) 

at dt 


is the absolute acceleration, and SI is the instantaneous speed of rotation. The gradients are 
taken relative to R. Taking 


P d V = p dv 

and equation 66 into account, equation 71 becomes: 


/ 


p [ V U (r) + u • W U (r) + V - w] dv 




T = 0 


(73) 


in which the terms of higher order are neglected. By applying Gauss’ theorem to equation 
73, the differential equation of tidal oscillations takes the form: 

pw = p(VU + u • VVU) + V • (A • T) (74) 

which contains the influence of the deviatory prestress. From equations 69 and 70, 

A•T = T + (l-^0)S-SXw-e•S (75) 
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and 


V • (A • T) = V • T + (1 + 0) V • S - (V • e - V 0) • S 


- (S • V) X w - (V • S) X w - e •• VS 


or taking into account 


V •e=V0-Vx<«) 


so that 


V*(A*T) = V*t+( 1 + 0)V*S + «XV‘S 
+ (VX«)*S-(S*V)X«-e**VS 

equation 74 becomes: 

pw = pV(U + tj/) + p u • W U + V* T+(l +0)V. • S 

(76) 

+ w X (V- • S) + (V X «) • S - (S • V) X - e • ♦ VS 


Note the existence of a coupling between the prestress and the vorticity. This coupling 
disappears if the prestress is reduced to hydrostatic pressure. Thus 

(k + N) (77) 

where is a constant and perturbations in the speed of rotation. Our primary 

interest is in diurnal and semidiurnal tides. Then N represents the nutation, astronomical 
and free. The typical term in N can be written in the form : 

N = c (i cos at - j sin at) 


thus, 


k ♦ N = 0 


(78) 
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In the system rotating with the Earth, the astronomical nutation is represented by a trigono- 
metrical series with the same arguments as in diurnal Earth tides (Reference 10). 

Substituting 

9R _ ^ 

3t ~ at 


into equation 72 and making use of equation 77 yields: 


a^u 3u dN 

w = —;r + 2S2 k X - S2 r X 

at^ " at ° dt 


(79) 


- V [0^ - (r X k) • (N X r)] - u . VV0^ 


where 


= — r • (I - k k) • r 


(80) 


is the centrifugal-force function. Assuming the static equilibrium in the “initial” state, 


where 


V • S + pVV = 0 


V = U + 0 


(81) 


(82) 


is the total gravitational force function. Using equations 78 and 81, differential equation 
76 of tidal oscillations becomes: 


P 



au \ 

+ k X ) = -p0VV + p(u • WV - w X VV) 

at f 


+ pV (0 - r * N k 


r) + V 


T + p 12 r X 

o 


dN 

dt 


(83) 


+ (V X w) • S - (S • V) X to - e • - VS 
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In the liquid core, the relation between the stress and strain tensors is assumed to be either 
Hooke’s law, 


T = X. 1 0 + 2u e = X 1 0 + (u V + V u) 

C 0 c C ^ ■' 

or, more generally, the Kelvin-Voigt law. 


where 


T = X I 0 + 2jLte = XI0 + n(uV + Vu) 


d 


b 

at 


(84) 


(85) 


( 86 ) 


For the oscillations of a ^ven frequency, o. 


X = Xg - iffXy.M = Me - 

Generally speaking, \ , and are functions of r. The values of X^ and n^aie supplied 

by Earth’s models. Knowledge of the viscosity of the liquid core is incomplete (Reference 
22). Estimates of the kinematic viscosity coefficient of the liquid core vary between lO'^m^s"^ 
and 10^ m^s'^ (References 21 and 23). In its main characteristics, the motion in the interior 
of the liquid core follows the laws of classical hydrodynamics (Reference 21). For this reason 
and to obtain at least preliminary information on the influence of the viscosity on the tidal 
oscillations, the numerical integration can be performed under the Stokesian approximation: 


p 


= V 


\ 

P 



V 


where is a constant. 

Equation 83 is the basic. Any other equation of tidal oscillations that appears in the present 
exposition is either a consequence or a paraphrase of equation 83. Of special interest are 
transformations that bring equation 83 closer to a form associated with the hydrostatic 
equilibrium. For example, making use of the identity. 


yields: 


u . VVV = <oXVV + V(u*VV)-e«VV 


/ a"u 



V • t - P0 vv 


( 88 ) 
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+ p V (i// + u • VV - r » N k » r) 


(88 continued) 


dN 

+ pQ. r X - (p e • VV + e ••VS) 

° dt 


+ (V X w) • S - (S • V) X w 


or, taking the identity, 


-pflVV + pV(u • V V) = -V •(pu)VV + V(pu •VV) 
-uX(VpXVV) 


into account yields: 




du 


, + 2S2 k X = V • T - V • (p u) VV 

° dt/ 


+ p V (i// - r*Nk*r) + V(pu* VV) - u X (Vp X VV) 


(89) 


+ p r X — - (p e • V V + e • • V S) + (V X w) • S - (S • V) X <o 

Considerable simplification can be achieved if the prestress is reduced to hydrostatic pressure, 

S = - pi 

then, 

p VV - Vp = 0 


and, as a consequence, 

VpXVV = 0 VpXVV = 0 

i.e., in the initial state, the equipotential surfaces-the surfaces of equal pressure and equal 
density— coincide. It is assumed that, in the initial state, the equipotential surfaces also 
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coincide with the surfaces of equal compressibility and rigidity. As a consequence of the 
hydrostatic hypothesis, 


pe*VV+e* 
(V X «) • S - (S • 
and equations 88 and 89 become, respectively. 


a^u 

at^ 


au 

+ 2Q k X 

° at 


+. V()// + u*VV~ r • 


/ a^u au 


au \ 
at } - 


+ V(pu •VV) + pV(\l>-Q^r 


• VS = 0 


V) X « = 0 


= — V • T - 0 VV 
P 


N k • r) + n r X 

o 


dN 

dt 


V«T-V«(pu)VV 


• N k • r) + p 12^ r X 


dN 

dt 


(90) 


(91) 


Equation 91 represents a generalization, with the effects of nutation and geostrophic force 
added, of the equation used by Alterman et al. (Reference 7) for computing free oscillations 
of the nonrotating spherically symmetric Earth. Equation 90 can serve as the foundation of 
the Molodensky theory of tides (Reference 3). Equation 90 is transformed here to a different 
shape, more comparable to the classical Molodensky equation. As previously stated, because 
of hydrostatic equilibrium, the equipotential surfaces— the surfaces of equal pressure and of 
equal density— coincide. It was also assumed that the surfaces of equal X and p remain coin- 
cident with the equipotential surfaces. Consequently, 

X' m' p' 

VX = — VV.Vju = ~ VV, Vp = VV (92) 

g g g 

where X', ji' , and p' are the derivatives in the direction of gravity V V, and g = | VV |. Using 
equations 92 and 85 and after some easy transformations: 

1 Xp' „ /\e\ p 

— V • T = — - 0 V V + V { ) +— (V^ u + ve) (93) 

P gP \ P / P 
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(93 continued) 


+ — [gu' + (Vu) • VV] 
gP 

By substituting equation 93 into equation 90, the Molodensky equation (Reference 3) is 
obtained: 


9^11 9u ^ , 

— — + 2£2 k X = V(^ + u»VV + - r • N k • r) 

at* ° at p ° 


p evv + — (V* u + A 6) + — [g u' + (V u) • VV] 
p gp 

dN 

+ r X — — 

° dt 


(94) 


where 


(? = 1 - 


gP^ 


(95) 


is the core stability parameter (Reference 1 7). If the viscosity in the liquid core is neglected, 
equation 94 becomes: 


a*u 9u Xd - 

— — + 2J2 k X = V (i// + U • VV + - n* r • IM k • r) 

at* ° at p ° 


dN 

0 V V + n r X 

° dt 


( 96 ) 


Finally, if the liquid core is in neutral equilibrium (i.e., the Adams-Williamson relation, 
P = 0, is satisfied (Reference 16)), the equation takes a very simple form: 



If only one nutational term is selected, complete identity with Molodensky can be achieved 
(except for notations): 


N = « (i cosat - j sin at) 

From equation 98, 

dN 

= - ak X N 

dt 


and, after some easy vectorial transformation. 


^ U.1V 

- V r • N k • r) + r X = - V (S2^ - a) r • N k • r] - 2J2^a k r 


Then, for example, tidal equations 90, 91, and 94 become, respectively: 


3^u 3u 1 

+2J2^kX = — vt-0 V 

3t^ 3t p 

+ V + u • V V - - a) r • N k • r] - a k r • N 


and 


/ 3^u 3u\ 

p ( + 2n^ k X — = V • T - V • (pu) vv 

V 3t^ 3t / 

+ V (pu •VV)+pV[</'-S2j^ (fljj - a) r • N k • r] 

- in a k r • N 

O 


3^u 3u r Xd 

+2t2„kX = VU + U-VV+ n (n 

o O I O O 

dt^ L ^ 

- /30VV +— (V^U + V d) +— [gu'+ (vu) 
P gP 

- 20,^ a k r • N. 


a) r • N k • r 


VV] 


(98) 


(99) 


• N (100) 


(101) 


( 102 ) 


(103) 
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TRANSFORMATION OF NUTATIONAL TERM 

From 

N = e (+ i cos at - j sin a t) 

= + i sin 0 cos + j sin 0 sin </? + k cos B 

k = + e cos 0 - e. sin 0 
r u 

and taking 

sin 0 S/ = +e'*'^sin $ 

S+i = -3e+‘'<’ sin 0 cos 0 S'/ = +3e"‘‘^sin 0 cos 0 

into consideration, 

r • N k • r = + — e + e'**'* S'/ ) 

6 ^ ^ 


and, taking equations 15 and 16 into account. 


V (r • N k • r) 



1 . 1 \ p+iat 


+ 




The following is obtained in a similar manner: 

k N • r = e r (e^ sin 0 cos 0 - sin^ 6 ) cos (at + ^p) 
and, taking equation 105 into account, 

e sin0 cos 6 cos (at + <^) == - — A'^ 

' ^ ^ 6 ^ 6 ^ 


Making use of equation 60 and equation 8 , 

sin^ 6 cos(at + <^) = ^ i ^ e'*'^ 

+ ( - — b/ + — i C * ) 

\ 12 2 4 ^ / 




^-iat 


( 104 ) 


(105) 


(106) 


(107) 


(108) 


(109) 


( 110 ) 
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Substituting equations 109 and 110 into 108, 


I 


kN • r = 




s+l 



e+iat 


+ er 





( 111 ) 


The following representations of the nutational term is obtained from equations 107 and 
111: 


- ff) V (r * N k • r) - 212^^ a k N • r 
= +er ( + - n" A+i +-^l iO 

\3°^ 6°2 2 °W 

-er ( + - a;' +- b;i +- in ac-/) 

\3°2 6°2 2 °V 


( 112 ) 


By introducing the tesseral-force function: 

ni = + — en^ r2 (S+I - Sl^ e'*"*) (113) 

i ^ O i i 

which arises from the variation of latitude, the nutational term can be written in the form: 


-n^ (n„ - a) V (r • N k • r) - 2n^ ff k N • r = vn^ 



ern„ff(C+i 


6+*^^ + C;^e-‘‘^) (114) 


Thus, the nutation, forced or free, modifies the lunisolar tidal potential and initiates the 
toroidal oscillations. For the astronomical nutation, the force function, Hj, is similar in 
form to the corresponding term in the lunisolar tidal-force function. It is possible to fuse 
both force functions together. As a result, small terms added to W would be expected to 
modify the coefficients of the Earth’s tidal elastic response slightly. Each particular frequency 
has its own unique Love and Shida numbers. 
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EXPANSION OF THE CORIOLIS FORCE 


In this section, the Coriolis force is expanded into spheroidal and toroidal components. In 
the mantle, the tidal oscillations are predominantly spheroidal. In the liquid core, however, 
the Coriolis force and the ellipticity produce both toroidal and spheroidal tides and the 
numerous couplings between spheroidal and toroidal tidal constituents. 

By assuming that the tidal oscillations can be expanded into a sum of terms of the form: 

“n =U„(r,t)A„ + V„(r,t)B„ + il„(r,t)C„ (115) 

and by taking equations 6 1 through 63 into account, the following is obtained for the 
typical term in the expansion of the Coriolis force: 


du 


2n„ k X — = in < -i 

° at “ 


(n + 1) (n + m) 9T 


+ i 


n (2n + 1) at 


n-l 


. (n+ l)(n + m) ax„ 

• .. A 

at 


2n+ 1 


„ + m 1 

. I + (n+l) I 

n(2n+l) J 


n-l 


av„ im / au„ av \ 

- i m A I + / B 

at " n(n+ 1) \ at at / " 


(116) 


m 


n (n + 1) 


+ i 


(n+ l)(2n+ 1) at 


ax^ n(n-m+l) ax„ 

C + i A„., 

at " 2n + 1 at 


n (n - m + 1) 31^^ 




(n+ i)(2n+ 1) yat at y ) 


Under the influence of the Earth’s rotation, the spheroidal (toroidal) oscillations produce 
toroidal (spheroidal) terms in the expansion of the Coriolis force. There is also a mutual 
influence between the tides of the (n - 1)^* , , and (n + 1)^* degrees. 

Substituting first n + 1 , then n - 1 , for n in equation 116 and retaining only the terms of 
the n^* degree. 


du 

in k X — = in 
° at ° 



(n - 1) (n - m) 
2n- 1 


ax 


n-l 


at 


m 

at 


(117) 
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(117 continued) 


(n + 2) (n + m + 1) 

«« • 

2n+3 at 

r (n-l)(n-m) 31^, m 3U„ m 3V„ 

.+ i I + • — - • '■ 

L n(2n-l) 3t n(n+l) 3t n(n+l) 3t 

(n + 2) (n + m + 1) 31^,^^ 1 

+ • B 

(2n + 3)(n+l) 3t J " 


r n- m 


(n-l)(n-m) 3V^_j 

m 

ax 

n 

■f 1 + 

n (2n - 1) 

at 

n (2n - 1) dt 

n (n+ 1) 

at 

n + m+ 1 


(n + 2) (n + m + 1) 



(n+l)(2n + 3) 

at 

(n + 1) (2n + 3) 

at J " 

) 


This equation shows the manner in which the tides of the n^* degree are being influenced 
by the tides of (n - 1)^* and (n + 1)^* degrees through the Coriolis force (References 6, 
18, 24, and 25). There is no couphng, however, between the tides of different orders. In 
particular, for oscillations of the first degree and order. 


3u 

2J2„ k X — = 2£2„ 

o at o 


+ + 


. 

av} 

9 3T} 

.'V 

at 

5 at , 

j av 

1 

1 

_ j. • 

^ 9T}\ 

2 at 

T 

10 at / 

au‘ 

9 

avj\ 

"at” 

~ 10 

at / 


(118) 


For diurnal tides of the second degree. 


3u 

2f2„ k X — = 2J2„ 
o o 


/ 1 3t} av^ 

\ 3 at at 


16 3T‘\ 

T IT/ 


(119) 
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(119 continued) 


+ i + - 


1 axj 1 au* 1 av‘ i6 ai^' 


6 at 6 at 6 at 21 at 


/ 1 au| 1 av{ 1 ai^ 4 au^ 16 av^' 

\ 6 at 6 at ^ 6 at 21 at 21 at 


For diurnal tides of the third degree, 


au ( / 4 aT‘ av^ 25 axj 

in k X. — = in < i ( + 

° at °| V 5 at at 9 at 


+ i I + 


+ + 


ax* 

1 

au* 

1 

av* 

25 

aTj\ 

at 

12 

at 

12 

at 

+ — 
36 


4 

av* 

1 


5 

au' 

25 

15 


12 

at 

36 

at 

36 



( 120 ) 


For the most important semidiurnal effects, 


au \ / ^^2 20 

2n k X = in <i - 2 — — -) 

° at ° ) \ at 7 at / 2 


/ 1 au2 1 av2 20 \ 

+ i ( + bI 

V 3 at 3 at 21 at / 2 


( 121 ) 


/ 1 5 3 U 23 20 avf 

+ I + ) 02 

V 3 at 21 at 21 at ' 2 



and for n = 3, m = 2, 


au 

2J2« k X 

® at 


ioax2 

2J2„ { i + 2 

° ' 5 at at 3 at 


( 122 ) 
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(122 continued) 


+ i 


(• 


2_K 

15 at 


1 

6 at 


6 at 


6 at > 


L±^ 3 . ±!Zi' 

y 15 at ' 15 at 6 at ’ 6 at ' 6 at 


The influence of free nutation on semidiurnal tides is considerably smaller than on diurnal 
tides. In the liquid core, the Coriohs force produces numerous couplings between the semi- 
divirnal tides of different degrees. 

If tides of a given frequency a only are of interest, then by replacing 


U by U , V„ by , etc. 

m*' m m*' m 


in equation 1 1 8 we obtain 


It 


(n - 1) (n - m) 

2 i n a k X u„ = - 2n„ a <1+ T. , - mV_ 


2n- 1 


n-l 


(n 


+ 2) (n + m + 1) ^ 1 ^ 

2n + 3 J " 


n 


- l)(n-m) 


m 


(2n-l) 


n-l 


u - 


m 


n(n+l) " n(n+l) " 


(123) 


(n + 2) (n + m + 1) 

+ — T ^ , B„ 

(2n + 3)(n+l) " ^ " 

[ n - m (n - 1) (n - m) m 

+ U , - — V , + T 

n(2n-l) "■* n(2n-l) n(n+l) " 

n + m + 1 (n + 2) (n + m + 1) 

(n+ l)(2n + 3) (n+ l)(2n + 3) 

In particular, for diurnal and semidiurnal tides up to the third degree, as a paraphrase of 
equations 118 through 1 22, 
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( 125 ) 
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INFLUENCE OF ELLIPTICITY OF EQUIPOTENTIAL SURFACES 


(128) 


In the frame of the present theory, it is assumed that the interior equipotential surfaces are 
ellipsoids: 


r |l + e (a) Pj (cos 0)J = a 

3 , 1 

P. (cos 6 ) = — cos^ 6 - — 
2 '' ^ 2 2 


( 129 ) 


where a is the eUipsoid’s mean radius, and e is its ellipticity. In the expansions only the first 
power of e is retained. By assuming that U (r) is the force function of self-gravitation associ- 
ated with a properly selected (or assumed) spherical model of the Earth, the following is 
obtained (Reference 6) for the potential of the ellipsoidal model with the centrifugal part 
included: 


V = U(r) + 


1 

I 


- c Pj (cos 6 ) 


( 130 ) 


where 


c 


= + 


2 

3 


re(r)g(r) + 


1 

1 


o 


g = - 


dU 

dr 


( 131 ) 
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From equation 130, for the undisturbed force of self-gravitation, 


VV = - 



dc 1 

— P2(cos0)Je^ 


— e . sin 6 cos 6 
r ® 


(132) 


where 


7 = 



(133) 


In the further exposition, it is assumed that the stability parameter, jS, with sufficient accu- 
racy can be considered either a constant (the case of the uniform stability) or a function of 
r. Taking into account 




= X„ s 

n n 


where 


X 


n 


dU 2 n(n+l) 

— a + — u 

dr r " r 


yields: 


VV = fix 


dc 3c 

. 7 A„ + — e P. (x) S + — Q 

n'n jj. r 2 n ^ n 


(134) 


(135) 


(136) 


where 


X = COS 0 


(137) 


From equation 3, 




(X) = - 


3 (n - m + 1) (n - m + 2) 


(2n+ l)(2n+3) 


n+2 


(138) 
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(138 continued) 


-3m^ +n 

+ A 

(2n-l)(2n + 3) " 

3 (n + m - 1) (n + m) 

^ T * (2n-l)(2n + l) 

From equations 16, 17, and 7, and taking the normalization conditions into account, 

1 g . (n + m- l)(n + m) ^ 

t'(i' +l)a^ ^ ' (n-l)(2n-l)(2n+l) 


+ 


- 3m^ + n 

n(n+l)(2n-l)(2n+3) 


6 


u,n 


(139) 


(n - m + 1) (n - m + 2) ^ 

^ (n + 2)(2n + 3)(2n+ 1) 

and 


^ \ . im (n - m + 1) ^ 

»'(p+l)a (n+ l)(n + 2)(2n+ 1) 


im (n + m) 

+ — 6 


n (n - 1) (2n + 1) 

and, making use of equation 60, the following expansion is obtained; 


(140) 


Q = 

n 


(n + m- l)(n + m) 
(n-1) (2n-l) (2n+l) 


B ^ + 

n-2 


- 3m^ + n 


n(n+l)(2n-l)(2n+3) 


(n-m+ l)(n-m + 2) ^ im(n+m) 

^ (n + 2)(2n+3)(2n+ 1) ^ (n-l)n(2n + I) 

im (n - m + 1) 

+ ^ C 

(n+ l)(n + 2)(2n+ 1) 

By substituting equations 138 and 141 into equation 136, 


(141) 
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dc 





(n + m - 1) (n + m) 
(2n- l)(2n+l) 


- 3m^ + n 3 

+ A„ + — 

(2n-l)(2n + 3) ” 2 


(n - m 
(2n 


+ 1) (n - m + 2) ^ ^ 

+ l)(2n + 3) " 


3c r (n + 
+ a — i — 

r L (“-!) 


m - 1) (n + m) 




- 3m^ + n 


B_ 


l)(2n-l)(2n+l) n(n+ l)(2n- l)(2n + 3) " 

(n - m + 1) (n - m + 2) ^ 

^ (n+2)(2n + 3) (2n+ 1) 

iin(n + ni) _ im(n-m+l) 


(n - 1) n (2n + 1) 


^ (n+ l)(n+2)(2n+ 1) 


.]x„ 


(142) 


By substituting n-2, n-l,n, n+l,n+2fornin equation 142 and by retaining only the 
terms of the n^* degree, 


-00VV = + pX^yA^ 

dc r 3 (n-m- l)(n-m) 

dr L 2 "-2 (2n-3)(2n-l) 

- 3m^ + n 

+ X 

" (2n-l)(2n+3) 

(n+m+ l)(n + m + 2) ~j ^ (143) 

+ 2 ^"+2 (2n+3)(2n+5) J " 

3c r (n - m - 1) (n - m) 

^ r L n(2n- l)(2n-3) 

n2 - 3m2 + n 

+ X 

" n(n+ l)(2n- l)(2n + 3) 
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^ (n + m + 1 ) (n + m + 2 ) ^ 

(n+ l)( 2 n + 3 )( 2 n + 5 )J " 


(143 continued) 


3 cm r n - m n + m + 1 ~| 

~~r L n(n+ l)( 2 n-l) ^ n (n + 1 ) ( 2 n + 3 ) J " 

In particular, for diurnal and semidiurnal tides up to the third degree, 

n= l,m= 1 , 

- ^0 VV = + PyX\ A} + 0 - Xj + — X‘ j A| 

3 | 3 c / 1 , 6 A , 9 c , , 

r V 10 ^ 35 V ^ 10 r 2 ^ 


(144) 


n = 2 , m = 1 


11 dc / 1 , 10 ,\ , 

-^evv = + /?7X‘Ai yX^+ 

c / 1 1 20 A , c / 1 , 2 A , 


(145) 


n = 3 , m = 1 


11 dc / 1 , 1 , 5 A , 

-P0VV = yPXlAl 7 + - X^ + - X‘j A^ 

p — (+ — X? + — X^ - — xA B* + — (+ X^ + X^) C 

^ r \ 15 ^ 20 3 22 V ^ r \ 10 ^ 36 


(146) 


n = 2 , m = 2 


- fidVY = + pyXlAl + p 


dc / 2 , 5 ,\ 

di V 7 2 7 


''a 


(147) 


c / 1 , 10 o\ o c 5 , . 

- 0 — {+ — Xl+ — X^ ) b 2 + ijS — . — X^ C; 

r \72 21 '’/ 2 r 7 ^ 


(148) 
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n = 3, m = 2 


-j30AV = + /37X2A2+/J 


1_ 

11 






(149) 


+ _ 


r 





+ 



c 


2 

3 


Because the spheriodal terms in equation 149 as produced by the ellipticity are extremely 
small, only the toroidal perturbative part might be significant. 

From equations 115 and 132, and taking equations 15 through 17 into consideration, 


• V V = - tU - — U P, (x)P;j' 

n n ' n n n 2 ^ ' n ^ 


t 3c 
— V 
r - 


^ dP"* (x) 

x(l - x^) 

dx 


+ m T 


„xP-(x)J 




(150) 


By taking equations 1, 2, and 7 into account, some easy transformations yield: 


ij = - yU S 

'n ' n n 


(n + m - 1) (n + m) r 3 dc 3c 1 

2 1- + U + (n + 1) — V Is . 

(2n-l)(2n+l) [_ 2 dr " t "J "'2 


3c m (n + m) 


T S 


2n+l " 


- 3m^ 


(2n-l) (2n 


+ n / dc 3c \ 

( + — U + — VIS 

n + 3) V dr " r "/ " 


(151) 


3c m (n - m + 1) 


T. S 


2n+l " 


(n - m + 1) (n - m + 2) / 3 dc 3c \ 

(2n+l)(2n+3) V T dT “ T ^"+2 
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Replacing n in equation 1 5 1 by n + 2, n + 1 , n, n - 1 , n - 2, and keeping only the essential 
terms (i.e., the terms of the degree n) yields; 

T? = u • VV = - yU S 

' ' n n 


( dc 


_3_ ,(n + m + 1) (n + m + 2) 

2 ’ (2n + 3) (2n + 5) 


U. 


n+2 


- 3m^ + n 
(2n-l)(2n + 3) 


(n - m - 1) (n - m) 
(2n - 3) (2n - 1) 


U „.2 


3c r (n + m + 1) (n + m + 2) (n + 
r L (2n+3) (2n + 5) 


3) 


V 


n+2 


m(n + m+ 1) 

+ — T 

2n+3 


- 3m^ + n 

^-l)(2n + 3) " 


m (n - m) 

+ — ^ 1 
2n - 1 


(n - m - 1) (n - m) (n - 2) 


n-2 


(2n-3)(2n- 1) 

In particular, for diurnal and semidiurnal tides up to the third degree, 


for n = 1, m = 1 , 


v[ = 


[I (* 


( 152 ) 


( 153 ) 
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(153 continued) 



144 , 9 , 3 , 

V* + — Ti V} 

35 ^ 5 2 5 ^ 


)] 


S 


1 

1 


for n = 2, m = 1 , 



(154) 


for n = 3, m = 1, 



(155) 


and for n = 2, m = 2, 


ril = - 



(156) 


+ 






s 


2 

2 


Let X, n, and p be the values of the elastic parameters and density for an ellipsoidal model 
of the Earth, and let , and p^ be their values for the corresponding spherical model. 

If the viscosity of the core is neglected, then X^^ , , and are functions of r. With the 

viscosity considered, they become the differential operators with respect to time: 


^ ® W 

2 

M = Mo + jM;re(r)P2(x) 


(157) 
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(157 continued) 


P 


Po +7p;re(r)P2<x) 


where 



0 

1 

1 

0 

dr ’ " 

° ■ dr ’ ' 

■ dr ’ 


From equations 157, 134, and 135, 



where 



(158) 


(159) 


Taking equation 3 into consideration, 

X X^ r (n - m + 1) (n - m + 2) 

— P=— -XS +S + — ^n+2 

P Pq " L (2n+l)(2n+3) 

2 - 3m^ + n (n + m - 1) (n + m) 

3 (2n - 1) (2n +3) " (2n - 1) (2n + 1) 



(160) 


Replacing n in the last equation by n - 2, n, n + 2 and keeping only the terms of the n^* 
degree. 


— 0 = — X S + s 
P P " " 


H 


2 - 3m^ + n 

3 (2n-l)(2n+3) 


(n - m - 1) (n - m) 

+ — — X 

(2n-3)(2n-l) 

(n + m + 1) (n + m + 2) 


n-2 


X + 

" (2n + 3)(2n+5) 


‘.J s. 


(161) 
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mill I II 


I I III I II 


llll■lll■H I mini III I I 


In particular, from equation 161 for diurnal and semidiurnal tides up to the third degree, 


n = 1, m = 1, 


X 

P 


12 




(162) 


n = 2, m = 1, 


— e — x‘ + s 
P Po 


( 2 20 \ 

+ — xl + -— xn s‘ 

21 ^ 63 '‘7 2 


(163) 


n = 3, m = 1, 


— e = — x‘ si + s 

P Po 


( '2 2 10 \ 

+ — Xi + — X* + — XU 

^ 15 ^ 15 3 33 sy 


(164) 


X 

P 


n = 2, m = 2, 
/ 4 


10 


— 8 = — X^ + s — X\ + — Xlj Si 

p 2 2 \ 21 2 21 2 


(165) 


n = 3, m = 2, 

X X„ 14 

— e= — x\sl +— sxlsl 


33 


■^5 3 


(166) 


SCALAR FORM OF DIFFERENTIAL EQUATIONS OF TIDAL OSCILLATIONS 

Substituting the expansions given in the three previous sections into the differential equation 
of tidal oscillations (e.g., into equation 103) and setting/Lx = 0 for the liquid core, for the 
general core of tidal oscillations of degree n: 


a^u 


at^ 


[• 


n . I (n-l)(n-m) 91^^ 

+ 212^ i 1+ - m 


2n- 1 


at 


av 

n 

at 


(n+2)(n + m+l) 

2n+ 3 at J 


(167) 


dc 

— + 
dr 


3 (n - m - 1) (n - m) 

= +07X +0— |+— -X 


2 (2n - 3) (2n - 1) 


n-2 
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- 3m^ + n 3 (n + m + l)(n + m + 2) 

+ — X + X + 

(2n-l)(2n+3) " 2 (2n + 3)(2n+5) dx 


n 

+ 2J2 

de 

(n-l)(n-m) dT^j (n + 2) (n + m + 1) 

+ + 

n (2n - 1) 9t (2n + 3) (n + 1) 

3c I (n - m - 1) (n - m) - 3m^ + n 

O I + i ^ 

r I n (2n - 1) (2n - 3) n (n + 1) (2n - 1) (2n + 3) 

(n + m+ l)(n+nH-2) ^ 1 ^ _1 
(n+ l)(2n + 3)(2n + 5) r 




where 


F„ = 


dc 


3 (n-m- l)(n-m) 

+ — L- 1 -^ 1 U , 

2 (2n - 3) (2n*- 1) 


n^-3m^ + n 3 (n + m + 1) (n + m + 2) 

+ U + — ^ U 

(2n-l)(2n+3) ^ 2 (2n+3)(2n+5) 


n+:^ 


4 [- 


(n - m - 1) (n - m) (n - 2) - 3m^ + n 

(2n-3)(2n-I) (2n - I) (2n + 3) 


(n + m + 1) (n + m + 2) (n + 3) 
(2n+3)(2n+5) 


] 


-t[ 


m (n - m) m (n + m + 1) 

2n- 1 2n + 3 


.] 


(167 

continued) 


( 168 ) 


( 169 ) 


40 



\ r (n - m - 1) (n - m) 

* * r (2n-3)(2„-I) ’‘-3 continued) 

,2 - 3m^ + n ^ (n + m+l)(n + m + 2)^ 1 

3 (2n - 1) (2n + 3) " ^ (2n + 3) (2n + 5) 

The nutational term 

ni = + - eS22 
2 6 ° 

must be added to the right-hand side of equation 169 for n = 2, m = 1, 

32 t^ r n-m 3U„.j (n-l)(n-m) 3V„.j 

20^ i + 

3t^ |_ n(2n-l) 3t n(2n-l) 3t 

n + m+1 3U„^i (n + 2)(n + m+l) 3V„^j m 3T„1 

(n + 1) (2n + 3) 3t (n+l)(2n + 3) 3t n(n+l) 3t I (170) 


3c 

+ /? — m 

I 


+ 


n-m 

n(n+ 1) (2n- 1) 


X 


n-l 


n+m+1 

+ — 

n (n + 1) (2n + 3) 



The nutational term, 


+ — 8 r 
2 ° 


must be added to the right-hand side of equation 1 70 if n = 1 . If interest is in tides of a 
given frequency only, taking equation 123 into account yields: 


- a- 


U -2 n a 

n o 


^ (n- l)(n-m) ^ 
2n - 1 


- m 


(n + 2) (n + m + 1) 
2n+3 



(171) 


dc 

= +/3yX„+/3 — 

dr 


3 (n - m - 1) (n - m) 

+ — — X 

2 (2n-3)(2n-l) ’ 
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(171 continued) 


- 3m^ + n 3 

+ X„ + — 

(2n-l)(2n + 3) " 2 


(n + m + 1) (n + m + 2) 
(2n+3)(2n+5) 


3r 


- V - 2 a 

n o 


[- 


m 


(U + V„) 
n(n+l) " " 


^ (n - 1) (n - m) ^ ^ (n + 2)(n + m+ 1) 

n(2n-l) (2n + 3)(n + l) 



3c r (n - m - 1) (n - m) - 3m^ + n 

r L n(2n- l)(2n-3) n („ + 1) (2n- 1) (2n + 3) " 

(172) 

(n + m+ l)(n + m + 2) ^ "j 

~ (n+ l)(2n + 3)(2n+5) ^ n 


- 0^ T + 2 a 

n o 


+ 


n- m 
n (2n - 1) 


U 


n-l 


(n- l)(n-m) ^ 
n(2n-l) 


n+ m + 1 
(n + 1) (2n + 3) 



(n + 2) (n + m + 1) m 

~ (n+ l)(2n.+ 3) * n(n+l) 

n-m n + m+1 

X , + 

n(n+l)(2n-l) n(n+l)(2n + 3) 



(173) 


Equations 167 through 173 show that there is a strong mutual interaction between U„> 
, and . The corresponding differential equations constitute a separate group. 
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Besides V^j, and j, some other terms enter into these differential equations, 
but they represent the effects of higher order. 


For example, for diurnal tides, 

/ 1 9T} dY\ 16 6Ti\ 

+ 2 ^2^ i 1+ " - — ■ 1 

9t^ V 3 9t 9t 7 at / 


= +PyX\ + P 


dc / 1 , 10 

— I + — xi + — ; 

dr \ 7 ^ 21 


9 Fi 

9r 


(174) 


a^V‘ / 1 9Tj 

+ 2 i (+ — 

at^ \ 6 at 


1 auj 1 ayi i6 ai^ 

6 at 6 at 21 at 


= +P 


r V 14 




F 


1 

2 


(175) 


a 2 


at^ 


/ 1 ai{ 

2 J 2 „i + 

\ 2 at 


3 au^ 9 a v‘ 

10 at 10 at 


(176) 


9 c , 1 .. . 

= + — P — xi + — e r f2 a 

10 r 2 2 ° 


a^ Tj / 2 a 

2f2^i (+ 

at^ \ 15 9t 


4 avj 1 3T‘ 

15 at 12 at 


5 au‘ 

36 at 


25 avA 

36 at / (177) 


= +/3— f+— + — xA 

r \ 10 2 36 




dc 

dr 


/ 1 , 10 
+ — U‘ + — 
V 7 2 21 



(178) 
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(178 continued) 


where 


c / , 3 , 12 , 100 

^ + s (+ — xj + — X* ) + 

„ * \ 21 2 63 4 / 2 


n‘ = +-^eJ2^ 


and 




+ 2f2_i 


dt* 


'(•7 


3T^ 3V^ 25 3T*'' 


3t 3t 9 3t 


= +Pyxl 


dc / 1 1 5 \ 3F^ 

+ ^ — (+ — X} + — xi + — xi ) + 

dr \ 5 * 5 3 11 7 3r 


3^V^ 

3t^ 


4 

3Tj 

1 

3U^ 

1 

3V» 

25 

15 

3t 

12 

3t 

12 

3t 

36 


c / 2 , 1 , 5 ,\ 1 , 

= +0— (+ X? + — X‘ X|) +— F* 

r \ 15 ^ 20 3 22 V r ^ 


d^T\ ^ 1 3T^ 1 3U} 1 3V} 4 3U‘ 16 3V|' 


3t^ 


6 3t 6 3t 6 3t 21 3t 21 3t 


= f 7 ’‘i 


(179) 


(180) 


(181) 
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q2j1 

at* 


- 2 S2„ i 

O 


/ 1 ax‘ 3 au* 9 av‘ 6 au* 36 av‘\ 

\ 2Q at 28 at 28 at ss at ss at / 

3c / 3 , 3 

= +/3 — (+ X* + — xM 

r \ 140 * 110 V 


(182) 


^ i dc / 1 - 1 , 1 

F. .yVl-- U _u|, ^U: 


5 

11 ® 


c/2, 3, 60, 6, 5. 

— V? + — Vi + — Vi + _ Ti + — Ti 

5‘ 5* 11® 5* 3'* 


(183) 


+ Xi + s (+ — X] + — Xi + — Xi 


15 


15 


33 


and, equations 121 and 122, 147 and 148 yield for main semidiurnal tides. 


/ av* 20 ax|\ 

\ at 7 at / 

, dc / 2 , 5 -\ apl 


a*u* / a V* 

+ 2 n i I - 2 

at* 


(184) 


a*v* 


at* 


, / I dvl 1 av* 

+ 2 ^2^ i I + 

V 3 at 3 at 


20 ax*\ 

21 at / 


(185) 


c / 1 , 10 ,\ 1 

= (+— X* + X*) +— F 

r V 7 * 21 r 
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at^ 


/ 1 3 2 

2S2^i l+- 

\ 15 3t 15 


dVl 1 dTl 1 dUl 5 3V2> 


+ — 


at 6 at 6 at 6 at 


= +0 


— (+— +-^ 

r \ 10 2 3 4/ 


(186) 




/ 2 , 5 A c / 6 , 30 , 50 

(-y -^7 v^)-7 {-T^l ^7 J 


(187) 


+ — + s 

Po 


(4 


10 > 


If only terms of a given diurnal or semidiurnal frequency, a, are of interest, then replacing 


U byU byV e+‘"*,etc. 

n n ’ n n ’ 


in equations 174 through 187 yields: 


u; + 


° (-7 -"T 

, dc / 1 , 10 ,\ dFi 


(188) 


(- Vi + n a (- — l\ - — T\ +— U 2 ) 

t " 3 o y ''2 o \ 3 I 21 3 3 2/ 


= + 


1 c / , 40 ,\ 1 , 

— P—\Xl Xi ) + — Fj 

14 r V 2 9 4/ r 2 


(189) 


(-a" + n„a)T} + (- y -7 V») 


(190) 
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9 c , 1 

= + — /J — Xi + — cJ2or 

10 r 2 2 

(-o^ +— n or)Tl + S2„a (+ — Uj vn 

^ 6 ° ® ° V 15 2 15 2/ 


(190 continued) 


= +0 


- (+ — x\ +— xA 

r \ 10 2 36 '*/ 


- a 


n *2v;) 


, dc / 1 , 1 , 5 ,\ dF| 

= +^yxi + p — (+ — x; +— Xl +— xi) + 

3 dr \ 5 1 5 3 11 V dr 

(-a^ + — n aWl + (- — T* - — Ti +— uA 

\ 6 ° / 3 ° \ 15 2 18 ^ 6 3/ 


c / 2 , 1 , 5 1 

= +j3— + — x! + — Xi - — Xi +— F 
r V 15 ^ 20 3 22 ®/ r 


t; . n„» (.i U1 - i VI - if u; V.) 

(.a^ +1 a) Tl +n^a (+ — V\ V\ U> -— vA 

\ \0 ° f * ° \ 14 3 14 3 55 ® 55 ®/ 

9 c / 1 , 1 \ 

= + __ ,3 _ (+ _ Xi + — Xi 

10 r \ 14 ^ 11 5/ 


(191) 


(192) 


(193) 


(194) 


(195) 
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Hnd, for the main semidiurnal tide: 


-a^Ul + Tl + 4V2) 

dc / 2 5 \ dp2 

+T«o‘^) vl + Tl +j u*) 


(196) 


(197) 


(- 



i2 o 

O 


Tt + 


S2 o 

O 








(198) 


The lunisolar tides, Uj , Vj , do not exist, and it appears that the influence of tides , 

Uj, V^, U 4 , and V| in equations 188 through 198 is negligible. Of course, this assumption 
requires a full confirmation by numerical integration, but at least it can be used at the 
beginning of the computational process. Differential equations 1 88 through 198 are then 
split into three separate groups that can be solved independently. 


TRANSFORMATION OF THE POISSON EQUATION 

The tidal forces cause a small redistribution of matter inside the Earth. This results in a 
small increment ^ of the interior geopotential that satisfies the Poisson equation: 


y 2 ^ Q y . (p u) 


(199) 


By assxuning the expansion. 




gm 
’^n n 


(200) 
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equation 199 can be reduced to an equivalent set of ordinary differential equations for v//™ 
that must be solved with the corresponding set (equations 167, 168, and 170) for the tidal 
displacements and i//™ . 

By substituting 


P = Po (r) + 9(r) ?2 

2 dp 

q = — r — e(r), x = cos0 

3 di 


(157) 


and 


u = U„ A„ + V„ B„ + i T„ C„ 
n n n n n n 


(116) 


into V • (p u) and by taking equations 15 through 18 into account, some easy vectorial 
transformations result in: 


V ‘ 0>u) = 




+ 3- 


+ 3 — V„ X (1 - x^) 
r " 


/ 

dq 

+ 

U + 

\ 

dr " 

pm 

n 

(X) 

dP™ 

) 

n 

\ g+im¥> 

dx 

) 


+ qX J P,(x) 


By taking equations 1,3, and 7 into account, the foregoing equation becomes: 




+ P. X_^ 


V •(Pu) = \^U„+p„X„y S„ 


2^(n + m)(n->-m-J) rdq ^ ^2(n+l)-5-vl S„, 

2 (2n+l)(2n-l) [dr " ^ " r "J 


(201) 


m(n + m) q 

+ 3 — — T„ S„ , 

2n + l r " 


49 



- 3m^ + n 
(2n- i)(2n + 3) 


\dr 


Un + + 

n ^ n 



(201 continued) 


m (n - m + 1) q 
+ 3 2 . T S„, , 

2n+l r " 



(n - m + 1) (n - m + 2) / dq 
(2n + 1) (2n + 3) (17 


Un + qx„ 



S 


n+2 


In equation 201 , replacing n by n + 2, n + 1, n - 2, and keeping only the terms of the n^^ 
degree yields: 


V 



+ P 



+ 


3 

2 


(n + m + 2) (n + m + 1) 
(2n+5)(2n+3) 



U 


n + 2 


+ q \+2 2(n + 3) 


— V 

j. n+2 


+ 3 


rn (n + m + 1) q 
2n + 3 r 


T 


n+l 


- 3m^ + n 
(2n-l)(2n + 3) 


\dr 


U + q + 3 — 

n ^ n ^ 



( 202 ) 


+ 3 


m (n - m) 
2n - f~ 



2 (2n-3)(2n- 1) Ldr 


"d 


3 (n - m - 1) (n - m) dq q 

+ - • ^ ^ 1-^ + qX „.3 -2(n-2)^ V„., 


Substituting equations 201 and 202 into equation 199 yields the differential equations 
satisfied by \p^ : 


d^^P„ 2 n(n+l) 


dr^ r dr 


\p^ = +47tG 


le 




(203) 
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3 (n 
+ — • — 

2 (2n 


+ m + 2) (n + m + 1) fciq 
>n + 5) (2n + 3) [dr 


+ 2(n + 3) 


iv 1 

I "*’J 


m (n + m + 1) q 

+ 3 • ^ — T_,, + 


- Sm^ 


2n + 3 r (2n-l)(2n + 3) 


+ n / dq q \ 

— U + qX + 3— VJ 

In + 3) V dr " " r 7 


m (n - m) q 
+ 3 • — T 


2n - 1 r 


n -1 


(203 continued) 


3 (n - m - 1) (n - m) dq q 

"T • 0„-3)0n-l) t 


]l 


If the ellipticity of equipotential surfaces is neglected, the standard equation for a spherical 
model is: 


dr^ r dr 


n (n + 1) /dp^ ' 

«//„ = +47tG U + p„X„ 

r \dr , 


(204) 


Some interesting conclusions can be derived from equations 203 and 204. It is clear that, 
if the Earth is assumed to be spherically symmetric and originally in hydrostatic equilibrium, 
then toroidal oscillations do not produce any change in the geopotential. This is a well- 
known fact, confirmed by the gravimetric observations on the surface of the Earth. In the 
mantle, the influence of the ellipticity of equipotential surfaces on the oscillations is very 
small, and the tides are predominantly spheroidal. The introduction of Love and Shida 
parameters as functions of r takes care of these two facts in the mantle. 

However, in the liquid core, the ellipticity and the dependence of the density on the latitude 
(i.e., the lateral inhomogeneity) cause the toroidal terms to be present in equation 203, and 
the toroidal oscillations do influence the variations of the geopotential if they are sufficiently 
large. In general, with the introduction of lateral inhomogeneities into the Earth’s model, 
their influence on the tidal oscillations can be expected to increase. 

For the main diurnal tides, the equation for the variation of geopotential becomes: 


dr^ r dr 


'Pi 


= +47tG 





+ 




(205) 
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(205 continued) 


f 


i (^1 

1 \dr 


Ul +qXi + 3-^ 





10 /dq 

^ 21 \ dT 


Ul + qXi + 10-^ 



If the influence of and is neglected, then equations 205 and 188 through 191 consti- 
tute a separate system to be solved for U^, V^, TJ , Tj, and 

For the main semidiurnal tides, equation 203 becomes: 


d^vl'j 2 di//2 

dr^ r dr 


6 


i/Zj = +47 tG 






2 

1 


/dq 

Vdr 


+ 


qXl +3 


q 

r 



30 q 

T 7 


t>2 

^3 


7 Vdr 




+ qx? + 



(206) 


Neglecting the influence of and results in a complete system consisting of equations 
206 and 196 through 198 to be solved for tides U^, V|, and Tj and for i//|. 

GENERALIZATION OF SOME MOLODENSKY EQUATIONS IN THE 
RESONANCE THEORY 

The close proximity of the period of the diurnal wobble of the liquid core to the sidereal 
day produces a resonance effect that influences the amplitudes of nutation and the value 
of the diurnal Love number. 

Molodensky and Kramer (Reference 3) have developed the theory of the resonance effect, 
assuming the neutral equilibrium of the liquid core and the validity of the Adams-Williamson 
condition (Reference 16). 

Several works recently appeared (References 13, 17, and 18) that consider the dynamical 
consequences of the departure of the liquid core from the neutral stability. In this section 
the Molodensky and Kramer (Reference 3) theory is extended and the model of the liquid 
core is considered, which departs from the neutral stability and, thus, does not satisfy the 
condition of Adams-Williamson (Reference 16). The theory is of the first order (i.e., the 
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influence of the ellipticity on the spheroidal tidal oscillations can be neglected), but the 
eccentricity shall be retained as a factor multiplying large toroidal oscillations and, in partic- 
ular, in the boundary conditions on the core-mantle interface. 

Molodensky and Kramer differential equations (Reference 3) are amended by perturbative 
terms proportional to the Pekeris and Accad index of stability of the outer core (Reference 
17). 

If a ~ then the toroidal oscillation, Tj , becomes large and, because of the small factor, 
- a, equation 190 is no longer suitable for determining tidal amplitudes. Equations 188 
through 191 must be replaced by a different combination of equations. Multiplying 
equation 190 by -1/3 and adding the result to equation 189 yields: 


1 1 2 
+ — a (a - 2S2 ) r T, + — e ai = F, + <l> 

o ^ 1 ^ o 2 


(207) 


where 




32 

8 / , 

14 

— 12 oT 

n a U, + - 

— CT 

21 o 3 

15 ° \ 

15 ° , 

8 

o c , 20 c 


35 

/3-X - — p—X. 
r ^ 63 r ^ 



(208) 


Equation 188 can now be rewritten in the form: 


2 dF2 

— n ctT, + = 

3 ° ‘ dr 


(209) 


where 

, / 32 \ dc / 1 10 \ 

^ = -a^U, + l+—T^^2V^j.pyX, - IS- ^ (210) 

For brevity and because the ambiguity does not arise, the superscripts and the factor, , are 
omitted in the exposition. 
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Eliminating F 2 from equation 209 by means of equation 207 yields: 


a(a - 2 £2^) 


d 


(rT,) + 2fi ctT, + e£2 ar = 3 

1^ o 1 o 



( 211 ) 


In this case, the following estimates of orders of magnitude are obtained: 


-o c 

O 


and, from equation 190, 


Cl -a 

O 


U V T 


eT, 


( 212 ) 


(213) 


and, from equation 135, 


X 


2 



(214) 


Consequently, 


1 d$ 

— - eCllJ^ (215) 

These estimates show that the terms on the right-hand side of equation 21 1 are less significant 
than the terms on the left-hand side. Furthermore, the form of equation 211 suggests the 
estimate, 


T 


1 


a 

- — r ~ e Tj 
2 ^ 


(216) 
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where a is a large constant (Molodensky and Kramer resonant parameter (Reference 3)). 
From equation 211, 


a 


o 


12 -a 

o 


(217) 


Setting 


- a 

2;; = e--^ a (218) 


and keeping only the essential terms in equation 1 54 yields: 

“2 --7 ♦7T,) (219) 

and the estimate 

77^ ~ egTj (220) 

From equation 178 the following can be derived with the same accuracy as that given 
previously: 


F, = 'P. 





( 221 ) 


and, from equation 190 




1 

2 S2 a 

O 



( 222 ) 
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and taking into account equations 216 and 218, 


V 


2 



1 C 

- /? 

2 n^ai 


( 223 ) 


It follows from equations 212 and 214 that the last term on the right-hand side of equation 
223 is of the order of /? e , and it must be retained if the influence of the departure of the 
core from the neutral stability is considered. 

From equations 135 and 223. 


1 - 


30c 


S2 or 

O 


1 d ^ 

,4 dr ^ 


10 

T 


(224) 


In agreement with Molodensky and Kramer (Reference 3), setting 


- K 


yields from equation 22 1 , 


X, = - ^ + K) 

and, substituting 

p _ 1 1 dp 

\ 1-|3 gp dr 


into equation 226, 


^2 

^ 1 - |3 gp dr 


1 dp 

(V^ + K) 


(225) 


(226) 


(227) 


(228) 
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Tl"' 


Eliminating U 2 and X 2 from equation 223 yields: 


V2 +7t,) 




c 1 dp 

• — — (t?2 + K) 


1 - j3 2S2 ai gp dr 


Combining equations 228 and 219 yields: 


dp 1 dp ^ c ^ 

(1 - p)p X 2 +- U, =--(K --Tj) 
2 dr ^ g dr r ^ 


(229) 


and, from this equation and equations 219, 224, and 228, 


_d 

dr 





r^ dp c 

+ - -^:(K — Tj) 

g dr r ^ 


^ r^ dp 
g dr 


2<3 

1 - 0 


. JL (rj + K) 
n^r 


(230) 


which is one of two basic Molodensky equations (in a slightly modified form) with the 
influence of the departure of the liquid core from the neutral stability included. 


Approximation 


Ti = 


1 

— ca 
2 


can be subtituted into equation 230. 

Neglecting the terms of the order e U 2 ~ in equation 205 yields: 


^^^^2 ^ 2 dip 2 
dr 2 r dr 


6 /dp 
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(231) 
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Substituting equation 219 into equation 231 for U 2 and into equation 228 for Xj, some 
simple transformations result in: 


^ 'I ' 2 2 ^'l'2 6 

dr2 r dr ,2 


^ 47tG dp /<^2 ^ 

g dr \ 1 - /3 



and, eliminating in favor of K by means of equation 225 yields: 


^2 r dr y g dr ^2 j 


+ 



47tG dp 
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P 

1 - 

- 


(t,2 + K) 


1 

-I 


(232) 


(233) 


in which the approximation, Tj = a r/2, can be used. To simplify equation 233, it is more 
convenient to introduce the new variable: 


Q = 


K 


— rT. 
3 ° ^ 


rather than K so that 
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+ 




(234) 
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From equations 207 and 218, 


-ia [(a-2J2„)rT, 




From equations 212, 220, 221, and 228, 


Q~ K ~ egTj ~ S2^rTj 


and, consequently. 


dr 2 r dr 



eQ 


r 


2 


and it can be neglected. Then, 


d^Q ^ 2 dQ / 47 tG dp ^ 6 \ 

dr" r dr \ g dr ,2 / ^236) 


4nG dp 
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g dr 



+ Q + 


— n^rT 'j 

3 ° 7 


or, in a slightly different form, 


d 

dr 


47tG 

= + 

g 


d 47 tG dp . 

r® T-(r'"Q) ^ 

dr J g dr 






(237) 


If the departure from the neutral stability is neglected, equations 236 and 237 become 
corresponding Molodensky equations. 

Equations 219, 228 through 230, 234, and 236 constitute a complete system for the res- 
onance case. It can be integrated only numerically. Solutions must be determined so as 
to satisfy the conditions on the boundary between the inner and outer cores and between 
the outer core and the mantle. 
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Because e satisfies Clairaut equation 


d^e 8jrGp de /SirGp 6\ 

— - + — + ( - — I e = 0 

dr^ g dr \ gr / 


(238) 


it can be shown that c satisfies the differential equation: 
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From equations 237 and 239, 
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(240) 


and combining equations 240 and 230 yields: 
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(241) 


which represents a generalization of Molodensky and Kramer equation 37 (Reference 3) to 
the case in which the liquid core deviates from the neutral stability. Integrating equation 
241 over the liquid core gives the extension of equation 39 of Molodensky and Kramer: 
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Equation 242 connects the resonance parameter, the ellipticity, and the index of stability 
with the boundary conditions. Euler’s equation for the angular momentum supplies an 
additional connection between the amplitude of nutation and the change in the tensor in 
inertia of the whole Earth, as produced by the small redistribution of mass (Reference 4). 
These two equations, together with the boundary conditions, constitute a complete set for 
determining necessary parameters of the theory, once the model of the Earth is given. 

CONCLUSIONS 

The present work provides the possibility of including the influence of the ellipticity and 
of the Coriolis force into the computations of tides in the Earth’s interior. It also deals with 
the mutual effects between the toroidal and spheroidal tides of different degrees. The 
Molodensky theory of resonance between the diurnal wobble and the diurnal astronomical 
tides of the liquid core is amended by including the possible deviation of the liquid core 
from the state of neutral stability. 

Plans for future work include the transformation of the boundary conditions and the con- 
struction of a model of tides in the outer core of the Earth. 

This model shall be tested against the variation of latitudes and positions of the pole and 
against the values of parameters of the Earth’s tidal elastic response (Love numbers) on the 
surface of the Earth. The analysis of the variation of latitudes and search for resonant fre- 
quencies was performed recently by Popov and Yatskiv (Reference 26) and Debarbat 
(References 27 and 28). It shall be repeated, making use of the observations of variation 
of latitudes and pole positions by modem technology such as laser ranging (Reference 11) 
and VLBI (Reference 29). The VLBI observations can also provide the values of Love 
number at each station and thus help to determine the inhomogeneities of the Earth’s tidal 
elastic response (Reference 29). 

Goddard Space Flight Center 
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